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We derive within a time-dependent scattering formalism expressions for both the current through 
ac-driven nanoscale conductors and its fluctuations. The results for the time-dependent current, its 
time average, and, above all, the driven shot noise properties assume an explicit and serviceable 
form by relating the propagator to a non-Hermitian Floquet theory. The driven noise cannot be 
expressed in terms of transmission probabilities. The results are valid for a driving of arbitrary 
strength and frequency. The connection with commonly known approximation schemes such as the 
Tien-Gordon approach or a high-frequency approximation is elucidated together with a discussion 
of the corresponding validity regimes. Within this formalism, we study the coherent suppression of 
current and noise caused by properly chosen electromagnetic fields. 
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I. INTRODUCTION 

The experimental success in the coherent coupling of 
quantum dotsiiS*^ has enabled measuring the transport 
properties of systems with a molecule-like level struc- 
ture. Recently, further progress in this direction has 
been attained by the reproducible measurement of cur- 
rents through molecules which are coupled to metal- 
lic leadsi^ Together with these experimental achieve- 
ments, new theoretical interest in the transport proper- 
ties of such nanoscale systems emerged. ^'^ One particu- 
lar field of interest is the interplay of the electron trans- 
port and excitations by an oscillating gate voltage, a mi- 
crowave field, or an infrared laser, respectively. Such ex- 
citations bear intriguing phenomena like photon-assisted 
tunnelingg-8-9 i0-"-i2.i3.i4.i5 ^nd the adiabatioiSiiLi2-La 

and non-adiabatio2Si2L22i2^ pumping of electrons. 

A prominent example for the control of quantum dy- 
namics is the so-called coherent destruction of tunneling, 
i.e., the suppression of the tunneling dynamics of a parti- 
cle in a double- well potential,^** in a two-level system, or 
in a superlattice.^^ Recently, coherent destruction of tun- 
neling has also been found for the dynamics of two inter- 
acting electrons in a double quantum doti2L2& Moreover, 
it has been demonstrated that a corresponding transport 
effect exists: If two leads are attached to the ends of a 
tunneling system, then a proper driving field can be used 
to suppress the current even in the presence of a large 
transport voltage. Moreover, in such a system the cor- 
responding shot noise level a priori can be controlled by 
proper ac fieldsi^ Within this work, we provide more 
details on this noise control scheme and also explore its 
limitations. 

An intuitive description of the electron transport 
through time-independent mesoscopic systems is pro- 
vided by the Landauer scattering formula'^^ and its vari- 
ous generalizations. Both the average current-^ and the 
transport noise characteristics'^!^ can be expressed in 
terms of the quantum transmission coefficients for the 
respective scattering channels. By contrast, the theory 
for driven quantum transport is less developed. Scatter- 



ing of a single particle by an arbitrary time-dependent 
potentials has been considered'^^''^^''^^ without relating 
the resulting transmissions to a current between elec- 
tron reservoirs. Such a relation is indeed non-trivial since 
the driving opens inelastic transport channels and, there- 
fore, in contrast to the static case, an ad hoc inclusion 
of the Pauli principle is no longer unique. This gave 
rise to a discussion about "Pauli blocking factors" 
In order to avoid such conflicts, one should start out 
from a many-particle description. In this spirit, within a 
Green function approach, a formal solution for the cur- 
rent through a time-dependent conductor has been pre- 
sented, e.g., in Refs. Is^ and|40 without taking advantage 
of the full Floquet theory for the wire. Nevertheless in 
some special cases like, e.g., for conductors consisting 
of a single level^i^ or for the scattering by a piecewise 
constant potential^iSS an explicit solution becomes fea- 
sible. Moreover, for large driving frequencies, the driv- 
ing can be treated within a self-consistent perturbation 
theory 

The spectral density of the current fluctuations has 
been derived for the low-frequency ac conduct ance^^*^ 
and the scattering by a slowly time-dependent 
potential^ For arbitrary driving frequencies, the 
noise has been characterized by its zero-frequency 
component A remarkable feature of the current noise 
in the presence of time-dependent fields is its dependence 
on the phase of the transmission amplitudes m^i^ By clear 
contrast, both the noise in the static case^^ and the cur- 
rent in the driven case^ depend solely on transmission 
probabilities. 

Within this work, we derive within a Floquet ap- 
proach explicit expressions for both the current and 
the noise properties of the electron transport through a 
driven nanoscale conductor under the influence of time- 
dependent forces. This generalizes recent approaches 
since the presented Floquet formalism is applicable to ar- 
bitrary periodically driven tight-binding systems and, in 
particular, is valid for arbitrary driving strength and, as 
well, extends beyond the adiabatic regime. The dynam- 
ics of the electrons is solved by integrating the Heisenberg 
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equations of motion for the electron creation/annihilation 
operators in terms of the single-particle propagator. For 
this propagator, in turn, we provide a solution within 
a generalized Floquct approach. Such a treatment is 
valid for effectively non-interacting electrons, i.e., when 
no strong correlations occur. Disregarding these inter- 
actions also implies that the displacement currents are 
not taken into account entirely. As a consequence, the ac 
component of the electrical current inside the nanocon- 
ductor may deviate from the particle current iSil^ 

This paper is organized as follows. After introducing 
in Sec. |n] a model for the leads and the conductor un- 
der the influence of external fields, we derive in Sec. IIIII 
for a situation with time-periodic but otherwise arbitrary 
driving general expressions for the current and its noise 
and establish a connection to a Floquet eigenvalue equa- 
tion. In Sec. IIVI we consider some special cases and 
approximations. Sectionals devoted to the influence of 
an electromagnetic dipole field on a conductor consist- 
ing of a few tight-binding levels. Situations with an ac 
transport voltage are addressed in Appendix ^ while in 
Appendix ^ we detail an alternative derivation which 
has been introduced in Ref. [s^- 



II. LEAD- WIRE MODEL 

We start out by introducing a model for the central 
conductor ("wire") under the influence of an external 
driving field like, e.g., a molecular wire subject to laser ra- 
diation or coupled quantum dots^''^ driven by microwaves 
or an oscillating gate voltage. The conductor is attached 
by tunneling couplings to external leads. The entire 
setup of our nanoscale system is described by the time- 
dependent Hamiltonian 

H{t) = iJwirc(^) + -ffloads + -^contacts, (1) 

where the different terms correspond to the wire, the 
leads, and the wire-lead couplings, respectively. We focus 
on the regime of coherent quantum transport where the 
main physics at work occurs on the wire itself. In doing 
so, we neglect other possible influences originating from 
driving induced hot electrons in the leads, dissipation on 
the wire and, as well, electron-electron interaction effects. 
Then, the wire Hamiltonian reads in a tight-binding ap- 
proximation with N orbitals \n) 

n.n' 

For a molecular wire, this constitutes the so-called Hiickel 
description where each site corresponds to one atom. The 
fermion operators c„, cjj annihilate and create, respec- 
tively, an electron in the orbital |n). The influence of 
an applied ac field with frequency fl — 2tt/T results 
in a periodic time-dependence of the wire Hamiltonian: 
Hnn'{t + T) = Hnn'{t)- The Icads are modeled by ideal 
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FIG. 1: Level structure of a nano-conductor with A'^ = 3 
orbitals. The end sites are coupled to two leads with chemical 
potentials /ii, and /i_R — /il 4- eV. 

electron gases, 

where c^^ (c^q) creates an electron in the state \Lq) 
{\Rq)) in the left (right) lead. The tunneling Hamilto- 
nian 

-^contacts = ^ {^Lqc\qC^ + ^flgcJ^qC^ 

+ h.c. (4) 

q 

establishes the contact between the sites \N) and the 
respective lead, as sketched in Fig. ^ This tunneling 
coupling is described by the spectral density 

T,{e)^2TTY,\Viq\''5{e~eq) (5) 

of lead £, £ = L,R. If the lead modes are dense, rg{e) 
becomes a smooth function. 

To fully specify the dynamics, we choose as an ini- 
tial condition for the left/right lead a grand-canonical 
electron ensemble at temperature T and electro-chemical 
potential Hl/R: respectively. Thus, the initial density 
matrix reads 

po oc e-(^'-^=-''^^^-''«^«)/^«^, (6) 

where Ni = '^tq'^tq is the number of electrons in lead I 
and fcsT denotes the Boltzmann constant times temper- 
ature. An applied voltage V maps to a chemical poten- 
tial difference /i^ — [i-L ~ with — e being the electron 
charge. Then, at initial time to, the only nontrivial ex- 
pectation values of the wire operators read (c|,g,C£g) — 

feieq)Sii'Sqq' whcrc fi{e) = (1 + exp[(e - ^£)/fc_BT])~i 
denotes the Fermi function. 

In our model Hamiltonian (Q, the leads are time- 
independent. Thus, it seemingly cannot describe ac 
transport voltages. Such a situation, however, can be 
mapped by a gauge transformation to one with time- 
independent chemical potentials as demonstrated in Ap- 
pendix ^ 
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III. SCATTERING APPROACH FOR 
TIME-DEPENDENT POTENTIALS 

Due to their experimental accessibility, the central 
quantities in a quantum transport problem are the sta- 
tionary current and the low-frequency part of its noise 
spectrum. Within a scattering picture of non- driven 
mesoscopic transport, both quantities can be expressed in 
terms of a transmission function T{E) which reflects the 
probability that an electron is transmitted from one lead 
to the other..?2. Due to energy conservation, the reversed 
process, occurs with equal probability. This is no longer 
true for driven systems and, consequently, the scatter- 
ing approach needs to be generalized. Thus, in this sec- 
tion, we derive expressions for the currents and its noise 
properties for the transport through the time-dependent 
system modeled above. In the so-called wide-band limit, 
the more compact derivation presented in Ref. ^3 be- 
comes possible, cf. Appendix IbI We will show that the 
average electrical current contains only transition prob- 
abilities and, thus, resembles a scattering formula. In 
clear contrast to the static two-terminal case, however, 
we will find that the noise depends in addition also on 
the phases of the scattering matrix. 



A. Charge, current, and their fluctuations 

To avoid the explicit appearance of commutators in the 
definition of correlation functions, we perform the deriva- 
tion of the central transport quantities in the Heisenberg 
picture. As a starting point we choose the operator 



We find below that for two-terminal devices Sg is inde- 
pendent of the contact £, i.e., Sl = Sr = S. 

The evaluation of the zero-frequency noise S directly 
from its definition 110(1 can be tedious due to the explicit 
appearance of both times, t and t — T. This inconvenience 
can be circumvented by employing the relation 



d 
dt 



{{Qm - {Qi{t)y 



dr Se{t, t - 



(11) 



which follows from the integral representation of Eqs. {Tj) 
and lO, Qt{t) — J^^dt' Ii{t'), in the limit to ~oo. 
By averaging Eq. Hll|) over the driving period and using 
S{t, t — t) = S{t — T,t), we obtain 



(|(AQ.^(t)) 



(12) 



where AQi = Qi ~ {Qt) denotes the charge fluctua- 
tion operator and (. . .)t the time average. The fact 
that the time average can be evaluated from the limit 
S — limtf,^-oo(A(5^(t))/(t — to) > allows to interpret 
the zero-frequency noise as the "charge diffusion coeffi- 
cient" . As a dimensionless measure for the relative noise 
strength, we employ the so-called Fano factoB^2*^ 



F = 



e|/| 



(13) 



QfXt) = eNi{t) - eNeito) 



(7) 



where / denotes the time-average of the current expec- 
tation value {Ii{t)). Note that in a two-terminal device, 
the absolute value of the average current is independent 
of the contact £. 



that describes the charge accumulated in lead i with re- 
spect to the initial state. Due to total charge conser- 
vation, Qi equals the net charge transmitted across the 
contact £; its time derivative defines the corresponding 
current 



(8) 



The current noise is described by the symmetrized cor- 
relation function 

Sdt,t')^i{[Am,AW)U) (9) 

of the current fluctuation operator AI({t) = Ig{t) — 
{Ii{t)), where the anticommutator [A, 5]+ — AB + BA 
ensures hermiticity. It can be shown that at long times, 
Sg{t, t') = Sg{t + T, t' + T) shares the time-periodicity of 
the driving.^^ Therefore, it is possible to characterize the 
noise level by the zero-frequency component of Si{t, t — r) 
averaged over the driving period. 



T 



dt 



dr Si{t, t — T). 



(10) 



B. Transition amplitudes 

In order to take the exclusion principle properly into 
account, we have formulated the transport problem under 
consideration in terms of second quantization. Neverthe- 
less, in the absence of interactions, both the current and 
its noise can be traced back to the solution of the cor- 
responding single-particle problem. Thus, our next step 
is to relate the expectation value and the variance of the 
charge operator to the transmission of electrons from 
one lead to the other. For that purpose, we start from 
the Heisenberg equations of motion 



(14) 



L/Rq^L/Rq^ 



(15) 

" = 2,...,iV-l. (16) 
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For these coupled linear equations, the formal solution 
ct'q'it) = Y.^tq'\U{t,t„)\tq)ctq{to) 

+ Y.'^i'q'\U{tM)\n) c„(io) 

n 

involves the propagator U{t,to) of the corresponding 
single-particle problem. We insert ((T7|) into {T)) and 
use the initial condition @ to obtain for the transfered 
charge at long times [i.e., in the limit to — > — oo, where 
all transients die out and, in particular, the second line 
in Eq. (|17|l becomes irrelevant] the expectation value 

(OlW) =e 5] {\{Lq'\U{t,to)\iq)f ~ SeLS,,,) Me,). 

9', 9/ 

(18) 

To symmetrize this expression, we eliminate the back- 
scattering terms, i.e., the contributions with £ = L, by 
employing the completeness relation 

q q n 

= Pl + Pb. + Pwire, (20) 

where Pl, Pr, and Pwire denote the projectors onto the 
states of the left lead, the right lead, and the wire, re- 
spectively. Then, from the time derivative of Eq. (|18|l . 
we find for the current through the left contact the result 

{Ihit)) ^e^{wLq'.,Rq{t)fBieq)~WRq,^Lq{t)fL{eq)} 

9:9' 

-eY,Wn,Lq{i)fL{eq) (21) 

n,q 

and mutatis mutandis for the current through the right 
contact. This expression already obeys the "scattering 
form" with the time- dependent transmission 

Ti,i{t, e) = 2^hY^Wi,q,^iq(t) 5{e - e,) (22) 

9,9' 

of electrons with energy e from lead t to lead At 
asymptotic times, the transitions from the lead state \lq) 
to the lead state \t'q') and the wire state \n) happen with 
the rates 

^/;r,'A(i) = lim ^\{eq'\U{tM)H\\ (23) 
Wn,iqii)^ lim ^\{n\U{tM)\^q)\^ ■ (24) 

The last term in the current (|21|) describes a periodic 
charging of the wire stemming from the external driving. 
With an average over one driving period, this contribu- 
tion vanishes and, thus, the dc current reads 

I = e^{wLq'MqfR{eq) - W R,q' ^Lqf L{eq)} , (25) 
9,9' 



with we'q'/q denoting the time average of the rate H23|l . 
Interchanging in Eq. L and R yields the negative 
current — /. Thus, as expected from total charge con- 
servation, the average current is, besides its sign, in- 
dependent of the contact at which it is evaluated. We 
emphasize that ()25|l obeys the form of the current for- 
mula obtained for a static conductor within a scattering 
formalism. In particular, consistent with Refs. J32 and 
Is^ no "Pauli blocking factors" (1 — fi) appear in our 
derivation. In contrast to a static situation, this is in 
the present context relevant since for a driven system 
generally WLq,B^' ^ WRq'.Lq, such that a contribution 
proportional to fL{eq')fR{eq) would not cance h^^i'^^ 

The zero-frequency noise S is conveniently derived 
from the charge fluctuation with the help of relation H12|l . 
Expressing the charge fluctuation by the Heisenberg op- 
erators ()17|) yields for the initial condition © after some 
algebra 

(AQi(t)) 

= E {Meq,)Meq)\{Rq'\U^PLU\Rq)\' 
q,q' 

+fL{eq')fR{eq)\{Lq'\U^PLU\Rq)\' 
+fL{eq')fL{eq) \{Lq'\UHPR^ + P^n-o)U\Lq)\' 
+fR{^q')h{eq) \{Rq'\U\PR + P^ire)U\Lq)\^}. 

(26) 

By using the completeness relation , we have achieved 
a form which is, besides the appearance of Pwire, sym- 
metric under exchanging L <-> i?. Here, C/ is a shorthand 
notation for U{t,to) and fi — I — fi- Taking the time 
derivative and averaging over the driving period yields 

S^e'J2{wkq,^RqfRieq')fRieq) 

q-.Q' 

+ Wtq',RqfL{eq')fR{eq) (27) 

+W^q'MfLi^q')fLi^q) 

+ W^q',LqfR(^q')fL{eq)}, 

where we have defined 

= ^Jmi^(^j^\{i'q'\UHt,to)Pe''Uit,to)\iq)f)^. 

The contributions in Eq. H26|) which contain the projec- 
tor Pwirc on the wire states do not contribute to the zero- 
frequency noise. This can be demonstrated by inserting 
for the propagator the explicit expressions H34|l and H35() 
which we derive in the next subsection. Interestingly 
enough, the noise S depends on both the diagonal and 
the off-diagonal elements of the projector WPinU. By 
contrast, the current depends only on the diagonal 
elements of this operator. As a consequence, in the pres- 
ence of driving it is not possible to express the noise solely 
by transmission probabilities; cf. Eq. (|45|l . below. 
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C. Lead elimination 

The evaluation of the rates w^'q'^iq and W^,^, in- 
volves the matrix elements of the time-evolution operator 
U{t,to) with the wire and the lead states. In the follow- 
ing, we eliminate the lead states and will find expressions 
for the rates that depend explicitly only on the propaga- 
tor for the wire electrons and the spectral density of the 
couphngs to the leads. 

We start from the Schrodinger equation for the prop- 
agator, ihdU{t,t')/dt = nit)U{t,t'), where n{t) is the 
single-particle Hamiltonian underlying Formal inte- 
gration with the initial condition U{t',t') — 1 results in 
the Dyson equation 

U{t, t') = Uo{t, t') - y* dt" Uo{t, i")Heo„tactsf/(i", t'), 

(29) 

where Uo denotes the propagator in the absence of the 
wire-lead coupling. We emphasize that due to the explicit 
time-dependence of the wire Hamiltonian, the integral in 
(EH is not a mere convolution. Using {i'q'\Uo{t,t')\£q) = 
Su'Sqqi eKp[~ieq{t — t')/h], we find for the transition ma- 
trix elements the relations 

{n\U{t,to)\iq) = - '-V;^ J' dt'e~i^^^''~'°\n\U{t,t')\m) 

(30) 

and 

{£'q'\U{t,to)\eq) 

= ,-(V*-,*o)|5^^,j^^, _ fdt' f'dt" 

^ Jta Jto 

xe*(v*'-^.*")(„,,|f/(i'^t")|„,)|^ (31) 

where ni denotes the wire site attached to lead £, i.e., 
Ti]^ = 1 and = N. 

At this stage, it is convenient to make use of the time- 
periodicity of the Hamiltonian, H(t) = 7Y(t + T). This 



has the consequence^ that U{t,t') = U{t + T,t' + T) 
and, thus, the retarded Green function 

G{t,e)^-- dTe^"U{t,t-T) = G{t + T,e) (32) 
" Jo 

can be decomposed into a Fourier series, G(i, e) = 
E^-oo e-''''"*G(''')(e), with the coefficients 

GW(e) / dte'^^^'G{t,e). (33) 
Jo 

Physically, G^*^' (e) describes the propagation of an elec- 
tron with initial energy e under the absorption (emis- 
sion) of \k\ photons for fc > (fc < 0). We emphasize 
that generally all sidebands k — — oo . . . oo contribute 
to the Green function H32(l and that, consequently, the 
fc-summations are unrestricted. 

After making use of Eqs. (|32|l and (|33|l . the transition 
amplitudes 1)30(1 and 1)3 l|l become 

(n|t/(t,io)|^g) = 1^4e-^^'(*'*°)5]e-''^-"*(n|GW(e,)|n,) 

fc 

(34) 

and 

{£'q'\U{t,to)\£q) 

^ k (35) 

respectively. Since below we restrict ourselves to asymp- 
totic times, to ~^ — cx), we have shifted the lower limit 
of the integrals accordingly. Moreover, in order to per- 
form the ^'-integration in Eq. H31|) , we have introduced a 
converging factor e''* 1^' and will finally consider the limit 
?7 ^ 0. 



I 

1. Average current 

For the further evaluation of the average current (|25|l . we insert the transition amplitude 135|l into H23|l . After 
taking the time derivative, averaging over time t, and considering the limit 77 ^ 0, we find 

fc 

and the corresponding expression for w^q' ,Lq- We have introduced the notation G„„' — (^n\G\n'). By use of the 
spectral density |SJ|, we replace the remaining sums over the lead states by energy integrals and obtain as our first 
main result the dc current 

^" = J E y - Tf^^^)^d^)] . (37) 

A;— — 00 
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where 

T^lI{^ =VL{e + khn)TR{e)\G'l},{t)\\ (38) 
T^nk^ ^TR(e + kh^l)TL{e)\Gf,{e)\\ (39) 

denote the transmission probabihties for electrons from the right lead, respectively from the left lead, with initial 
energy e and final energy e + kMl, i.e., the probability for an scattering event under the absorption (emission) of |fc| 
photons if A; > (A; < 0). 

For a static situation, the transmissions T''^^[e) and T^^{e) are identical and contributions with fc 7^ vanish. 
Thus, it is possible to write the current (|37|l as a product of a single transmission T(e) and the difference of the Fermi 
functions, /_R(e) — fhi^)- We emphasize that in the driven case this is no longer true. 



2. ac current 



Although below we focus on the computation of dc currents, we here continue the derivation of the transport 
quantities by presenting explicit expressions for the ac currents. We restrict ourselves to (/^(t)) since (/^(i)) simply 
follows by proper index replacements. Evaluating (/^(t)}, we consider also the last term in Eq. 121|l which describes a 
periodic charging/discharging of the wire. Apart from the time average we perform the same steps as in the derivation 
of the dc current and obtain 



(iLit)) = jjde {Tifl(t,e)/K(e) - T«i(i, e)/i(e)} - gi(i) 



(40) 



where 



n k 



(41) 



denotes the charge oscillating between the left lead and the wire. Obviously, since qL{t) is time-periodic and bounded, 
its time derivative cannot contribute to the average current. The corresponding charge arising from the right lead, 
qR{t), is a priori unrelated to qhit)] the actual charge on the wire reads qhit) -\- qnit). The time-dependent current 
is determined by the time-dependent transmission 



TLBit, e) =Tn{e) Re ^ e"^^^^* Gf^+'^(e) [Gfi,\e)\ * T^ie + k'hn) + -V de'- 



k,k' 



e' - e - k'hn 



(42) 



The corresponding expression for TjiL{t, e) follows from the replacement (L, 1) ^ (_R, N). Note that in the wide-band 
limit (e) = F^, £ = L, R, the contribution from the principal value integral vanishes. 



3. zero-frequency noise 

In order to obtain the zero- frequency noise S, we evaluate the rates W/,'^, This step is performed along the 
lines of reasoning for the evaluation of Wi'q'/q (although the actual calculation is far more tedious): We insert the 
transition amplitude (|35|l into H28|l . take the derivative with respect to t, and average over one driving period. Finally, 
we employ the relation lim^^o 4??[(e — a — iri){e — b + irj){e' — 6 — irj){e' — a + iri)]~^ = {2TT)^S{e — a) 5{(' — b) 5{a — b) 
to perform the limit 77 — > and find 



= |yi^J^^|5]r,.(e,-Kfc'ni^)[G('^;-*^H6,0]*G(^;^(£,) ^5{eq,-eq-khn), 



Rq',Rq 



(k'-k). 



(43) 



|Vfl,P5]|5]r4e,-t-A:'ftl7)[Glf-'=)(e,0]*Gl'^^(e,)-zGl'^(e,) ^ 5{eq, - tq - kh^l). (44) 
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The corresponding expressions for W^^, and Wj^^, follow from the replacement {L, 1) ^ {R,N). Inserting these 
into the noise expression H27|l we arrive at our central result 



fc k' 

+Tn{e,)TL{e)\Y,TL{e,,)Gf^-'\e,)[d^\e)Y 



iG 



i-k), 
IN 



iek)ffLie)fRiek)} 



(45) 



We have defined Ck 
density ©. 



+same terms with the replacement {L, 1) ^ {R, N). 
e + khfl and replaced the sums over the lead states by energy integrations using the spectral 

I 



D. Wide- band limit and Floquet theory 

In order to evaluate the expressions for / and S fur- 
ther, we derive an eigenfunction representation for the 
Green function. It is well-known that beyond the adia- 
batic limit, the eigenfunctions of the Hamiltonian are not 
of particular use — rather a proper basis is provided by 
a Floquet ansatzi^SiSi^ 

Let us start from the Schrodinger equation for the 
propagator. 



for n = 2 
d 



,7V - 1, and 



h- {ne \U{t, to)\n')^Yl i^"\U{t, to)\n') 



J2v;,{iq\Uit,to)\n'), (47) 



where rii is defined by = 1 and ur = N. To eliminate 
the lead states in the second line of Eq. H47|l. we insert 
(|Sn|l and replace by use of the spectral density @ the 
sum over the lead states by an energy integral. Then the 
last term in Eq. (|47|l becomes 



/der,(e) / dt'e-''^*'~'<^^/Hne\U{t,t')\n'). (48) 



Within the present context, we are mainly interested in 
the influence of the driving field on the conductor and not 
in the details of the coupling to the leads. Therefore, we 
choose for Ti{e) a rather generic form by assuming that in 
the relevant regime, it is practically energy-independent. 



(49) 



This so-called wide-band limit allows further progress 
since we now can perform in Eq. (|48|l the remaining en- 
ergy integration to obtain hS{t' — to) and, consequently. 



Eq. llT7|) becomes 
d 



n-{n,\U{tM)\n') =Y,Hn,n"{t){n"\U{tM)\n') 



-T,{n,\U{t,^)\n'). 



(50) 



Equations (|46() and (|50|l . together with the initial condi- 
tions {n\U{t,t)\n') = 6nn' fully determine the propaga- 
tor. Solving this linear set of equations is equivalent to 
computing a complete set of solutions for the equation 



^H-m)) 



where the self-energy 



+ \N)^{N\ 



(52) 



results from the coupling to the leads. 

Equation l|51|) is linear and possesses time-dependent, 
T-periodic coefficients. Thus, it is possible to construct 
a complete solution with the Floquet ansatz 

\tpa{t)) =exp[{-iea/fl- 'Ja)t]\Ua{t)), (53) 

oo 

|wa(i)) = \ua.k)e^pi-iknt). (54) 

A:— — oo 

The so-called Floquet states \ua{t)) obey the time- 
periodicity of TCwiieit) and have been decomposed into 
a Fourier series. In a Hilbert space that is extended by 
a periodic time coordinate, the so-called Sambe space;^ 
they obey the Floquet eigenvalue equationSiS 

(Hwivcit)-iT,~~ih-^^\ua{t)) = {ea-ih"/a)\ua(t)). (55) 

Due to the Brillouin zone structure of the Floquet 
spectrum^^Si^^i^ it is sufficient to compute all eigenvalues 
of the first Brillouin zone, —hil/2 < Eq < hVl/2. Since 
the operator on the l.h.s. of Eq. (|55l) is non-Hermitian, 
the eigenvalues Cq, — i^7Q are generally complex valued 
and the (right) eigenvectors are not mutually orthogo- 
nal. Thus, to determine the propagator, we need to 
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solve also the adjoint Floquet equation yielding again the 
same eigenvalues but providing the adjoint eigenvectors 
It can be shown that the Floquet states \ua{t)) 
together with the adjoint states \u'^{t)) form at equal 
times a complete bi-orthogonal basis: — 
Saf3 and J2a l'"ct(^))(^Q (01 ^1- A proof requires to ac- 
count for the time-periodicity of the Floquet states since 
the eigenvalue equation H55|) holds in a Hilbert space ex- 
tended by a periodic time coordinatei ^^'^^ 

Using the Floquet equation (|55|l . it is straightforward 
to show that with the help of the Floquet states \ua{t)) 
the propagator can be written as 

U{t,t')=J2e~'^'-^'-''-^^'-'"^\ua{t)){ut{t% (56) 

a 

where the sum runs over all Floquet states within one 
Brillouin zone. Consequently, the Fourier coefficients of 
the Green function [cf. Eq. H33|) ] read 



spin-dependent interactions, we have time-reversal sym- 
metry, WLq' ,R.q = WRqXqi and, sccoud, all sidebands with 
fc ^ vanish, i.e., r^2(e) = T^lr{^) = 4,oT(e), where 



T{e) = TL{e)Tn{e)\GM^)\' 



(59) 



and G(e) is the Green function in the undriven limit. 
Then the current assumes the known form 



(60) 



I F dt 



^ Jo ^ JO 

\Ua,k'+k){ul^k'\ 



Moreover in a static situation, the relation^SiS 

|ri(e)Gn(e)+*P = l-T(e), (61) 

allows to eliminate the backscattering terms in the second 
line of Eq. H45|l such that the zero-frequency noise can be 
expressed solely in terms of the transmission to read^ 



(62) 



E E 

a k' — — QO 



- 2 r 

'dTe'^^"'U{t,t- r) (57) ^o = ^ jde[T{e) [fL{e)h{^) + /i?(e)/fl,(e)] 

+T(e) [l-r(e)] [/^(e)-/i(6)f }. 



e — (e^, + k'hVL — ih'ya) 



(58) 



Inserting them into Eqs. %\7\ and (|45|l yields explicit ex- 
pressions for the current and the noise, respectively. 



IV. LIMITING CASES 

In the previous section, the dc current and the zero- 
frequency noise have been derived for a periodic but oth- 
erwise arbitrary driving. Within the wide-band limit, 
both quantities can be expressed in terms of the so- 
lutions of the Floquet equation (|55|) . i.e., the solution 
of a non-Hermitian eigenvalue problem in an extended 
Hilbert space. Thus, for large systems, the numerical 
computation of the Floquet states can be rather costly. 
Moreover, for finite temperatures, the energy integration 
in the expressions 137(1 and (|45|l has to be performed nu- 
merically. Therefore, approximation schemes which allow 
a more efficient computation are of much practical use. 

Before introducing various approximation schemes for 
the wire propagator, we discuss two particular cases for 
which current and noise assume more intuitive expres- 
sions. In doing so, we define quantities to which we will 
refer later in this section. 



A. Static conductor and adiabatic limit 

For consistency, the expressions H37f) and (|45|l for the 
dc current and the zero- frequency noise, respectively, 
must coincide in the undriven limit with the corre- 
sponding expressions of the time-independent scattering 
theory. This is indeed the case since the static situation 
is characterized by two relations: First, in the absence of 



For zero temperature, the terms in the first line vanish 
and pure shot noise remains. In contrast, for zero volt- 
age, fn = /l and the terms in the first line constitute 
equilibrium quantum noise. Obviously if both voltage 
and temperature are zero, not only the current but also 
the noise vanishes. In the presence of driving, this is no 
longer the case. This becomes particularly evident in the 
high-frequency limit studied in Sec. II V El 

It is known that in the adiabatic limit, i.e., for small 
driving frequencies, the numerical solution of the Flo- 
quet equation (|55|l becomes infeasible because a diverg- 
ing number of sidebands has to be taken into account. In 
more mathematical terms, Floquet theory has no proper 
limit as Oi^ The practical consequence of this is 

that for low driving frequencies, it is favorable to tackle 
the transport problem with a different strategy: If Ml 
is the smallest energy-scale of the Hamiltonian one 
computes for the "frozen" Hamiltonian at each instance 
of time the current and the noise from the static expres- 
sions H60f) and H62|) being followed up by time-averaging. 



B. Infinite voltage 

Many phenomena can be discussed in the limit of very 
large (practically infinite; subscript oo) voltages such 
that fR ~f 1 and ^ in the relevant energy range. 
Then, the dc current (|37|l becomes 

loo^^Y. j deTL{e + kh^l)TR{e)\Gf^{e)\\ (63) 



In the zero- frequency noise H45|l. only the contribution 

^q' ,q^^ Rq' ,Lq 



with /fl/i remains, thus, ^oo = e^Z^n'o^w. 
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^'^Y^q,q'i'^Rq',Lq " W^Lq',Lg)- derive this expression, 
we again have used the completeness relation (|19|l and 
the fact that terms containing the projector on the wire 
states do not contribute to time averages. Expressing 
WRq\Lq and W^^, by the Green functions yields 



5'cx) = e/oo--7- 2_, j deV L{tk)Y hie) 

k 



(64) 



where efc = e + khfl. These expressions make explicit 
that lac > and Sqo < sloe- Consequently, for infinite 
voltage the Fano factor H13() cannot exceed unity. 



C. Weak wire-lead coupling 

In the limit of a weak wire-lead coupling, i.e., for cou- 
pling constants Ti which are far lower than all other en- 
ergy scales of the wire Hamiltonian, it is possible to derive 
within a master equation approach a closed expression 
for the dc current The corresponding approximation 
within the present Floquet approach is based on treating 
the self-energy contribution — iE in the non-Hermitian 
Floquet equation (|55|l as a perturbation. Then, the ze- 
roth order of the Floquet equation 

(Wwhc(t) - ih^) \Mt)) = I0a(i)), (65) 

describes the driven wire in the absence of the leads, 
where |(/>a(i)} = J2k^^'Pi~^^^^)\'f'a,k) ai'e the "usual" 
Floquet states with quasienergies e°. In the absence of 
degeneracies the first order correction to the quasiener- 
gies is —ifvy]^ where 



1 



dt 

r 



(0„(t)|E|0,(i)) 



(66) 
(67) 



Since the first-order correction to the Floquet states 
will contribute to neither the current nor the noise, the 
zeroth-order contribution \ua{t)) = \u^{i)) = \(j)a{t)) is 
already sufhcient for the present purpose. Consequently, 
the transmission H38(l assumes the form 



E 

a,l3,k',k'' 



(jV|0a,fc')(0a,fc'+fc|l) 

E - (eo k'hn + ih-ii) 

{'^\4>fi,k"+k){4>p,k"\N) 

■ e - (fO -I- k"Ml - ih-yl) 



(68) 



and T^l(e) accordingly. The transmission (|68|l exhibits 
for small values of sharp peaks at energies e° -I- k'hn, 
and -I- k"hn with widths h'y^ and h'yl. Therefore, 



the relevant contributions to the sum come from terms 
for which the peaks of both factors coincide and, in the 
absence of degeneracies in the quasienergy spectrum, we 
keep only terms with 



a = k' = k" 



(69) 



Then, the fraction in (|68|l is a Lorentzian and can be 
approximated by 7r(5(e — e° — k'hfl) provided that 
is small. Consequently, the energy integration in H^Tfl can 
be performed even for finite temperature and we obtain 
for the dc current the expression 



E 



^LakVRak 



a,k,k 



'-[fRiel + k'hn)-fL{el + khn)]. 



The coefficients 

Tiafe =rL|(l|(/iQ_fc}p, 
^Rak =^TR\(N\(j)a,k}f, 



(70) 

rLo.=Y.rLo.k, (71) 

k 

r_RQ = ^ r^jctfe , (72) 



denote the overlap of the fcth sideband |0a,fc} of the Flo- 
quet state \(j>a(t)) with the first site and the last site of the 
wire, respectively. We have used 2?i7^ = TLa+^Ra which 
follows from (|67|l . Expression H70() has been derived in a 
prior workS within a rotating-wave approximation of a 
Floquet master equation approach. 

Within the same approximation, we expand the zero- 
frequency noise l|l5|l to lowest order in F^: After inserting 
the spectral representation (|58|l of the Green function, we 
again keep only terms with identical Floquet index a and 
identical sideband index k to obtain 



^_e- V- ^Rak'fRje^c+k'hn) ^ 2 p r (,o,,f,a\ 

a,k,k' (^La+i-Ra) 

+ {Tl„ + rj,jTLakfL{e°^+khn)} 

+ same terms with the replacement L ^ R. 

(73) 

Of particular interest for the comparison to the static 
situation is the limit of a large applied voltage such that 
practically /r, = 1 and /l = 0. Then, in Eqs. H7U|I and 
(|73|l , the sums over the sideband indices k can be carried 
out such that 

^LoTro 



=- y ■ 



(74) 



(75) 



These expressions resemble the corresponding expres- 
sions for the transport across a static double barrier.— 
If now Tlo = TflQ for all Floquet states \(j)a{t)), we find 
F — 1/2. This is in particular the case for systems obey- 
ing reflection symmetry^ In the presence of such sym- 
metries, however, the existence of exact crossings, i.e. de- 
generacies, limits the applicability of the weak-coupling 
approximation. 
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D. Homogeneous ac driving 

In many experimental situations, the driving field acts 
as a time-dependent gate voltage, i.e., it merely shifts 
all on-site energies of the wire uniformly. Thus, the wire 
Hamiltonian is of the form 



^wirc 



(76) 



where, without loss of generality, we restrict f(t) to pos- 
sess zero time-average. A particular case of such a ho- 
mogeneous driving is realized with a system that consists 
of only one level42iSiiS^ Then trivially, the time and the 
position dependence of the Floquet states factorizc and, 
therefore, the dc current can be obtained within the for- 
malism introduced by Tien and Gordon^ Here, we es- 
tablish the relation between such a treatment and the 
present Floquet approach. 

Since the time-dependent part of the Hamiltonian is 
proportional to the unity operator, the solution of the 
Floquet equation is, besides a phase factor, given by 
the eigenfunctions \a) of the static operator Hq — lE, 



Kit)) 



-iF{t)\ 



(77) 



where (Tig — iTi)\a) — (e^ — ifi'ya)\oi) and dF{t)/dt = 
f(t)/h. The quasienergies (e^ — itv^a) coincide with the 
eigenvalues of the static eigenvalue problem. Note that 
F{t) obeys the T-periodicity of the driving field since the 



J 



time-average of f{t) vanishes by definition. Thus, the 
phase factor in the Floquet states H77(l can be written as 
a Fourier series, 



-^F{t) 



E 



-ik^t 



(78) 



and, consequently we find \ua^k) = o,k\oi) and the ad- 
joint states accordingly. Then, the Green function 
becomes 



GW(e) = ^afe,+fea*, G{e-k'hn), 



(79) 



where G{e) denotes the Green function in the absence of 
the driving field. Inserting H79|l into (|37|l and employing 
the sum rule Y.k' <^k''^k'+k = 4,o, yields 



E 



h 



deT(e-fcm)[/fl(e)-/i(e)], (80) 



where T(e) is the transmission in the absence of the driv- 
ing. This expression allows the interpretation, that for 
homogeneous driving, the Floquet channels contribute 
independently to the current /. For the special case of a 
one-site conductor and a sinusoidal driving, this relation 
to the static situation has been discussed in Refs.l6lland 

El 

Addressing the noise properties, we obtain by inserting 
the Green function H79|l into H45|) the expression 



5-^E [de{\j2<^l'+kak'T{e-k' 

1. ■J 1.1 



m }R{e)fR{e + khn) 
+TlTr\ J2 al,+kak'GiN{e - k'hn)[TLGl^{e ~ k'm) - i\ |'/L(e)/fl(e + 

k' 

-I- same terms with the replacement (L, 1) ^ (i?, A^)|. 



(81) 



While the term in the first line contains only the static 
transmission at energies shifted by multiples of the pho- 
ton energies, the contribution in the second line cannot be 
brought into such a convenient form. The reason for this 
is that the sum over k' inhibits the application of the re- 
lation (|61|l . As a consequence, in clear contrast to the dc 
current, the zero- frequency noise cannot be interpreted 
in terms of independent Floquet channels. Only in the 
limit of large driving frequencies, we find below that the 
channels become effectively independent and reduces 
to an expression that depends only on the transmission 
in the absence of the driving and the Fourier coefficients 
afe, cf. next subsection. 

Expressions for the dc current and the noise that 
depend only on the static transmission have been de- 



r 



rived by Tucker and Feldman®^'^'* within a Tien- Gordon 
approach The central approximation of this approach 
is the description of a time-dependent chemical potential 
by an effective electron distribution. While this yields 
the correct expression (|Sn|) for the dc current, it does not 
capture the interference terms in the noise formula H81|l . 
This reveals that a Tien-Gordon-like approach yields the 
correct dc current while for the noise (and other higher- 
order correlation functions) it is only valid in a high- 
frequency limit. 

For large voltages where /l = and fR = l, the sums 
over the Fourier coefficients in Eqs. H80() and H81I) can be 
evaluated with the help of the sum rule ^j,, a'^,ak'+k = 
Sk^o- Then both the dc current and the zero- frequency 
noise become identical to their value in the absence of the 
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driving. This means that for a sufficiently large transport 
voltage, a time-dependent gate voltage has no influence 
on the average current and the zero- frequency noise. 

E. High-frequency driving 

Many effects occurring in driven quantum systems, 
such as coherent destruction of tunneling24 or current 
and noise controlj^SiSi are most pronounced for large ex- 
citation frequencies. Thus, it is particularly interesting 
to derive for the present Floquet approach an expansion 
in terms of l/f2. Thereby, the driven system will be ap- 
proximated by a static system with renormalized param- 
eters. Such a perturbation scheme has been developed 
for two-level systems in Ref. [s^ and applied to driven 
tunneling in bistable systemsS^ and superlatticesiS& For 
open quantum system, the coupling to the external de- 
grees of freedom (e.g., the leads or a heat bath) bears ad- 
ditional complications that have been solved heuristically 
in Ref. 44 by replacing the Fermi functions by effective 
electron distributions. In the following, we present a rig- 
orous derivation of this approach based on a perturbation 
theory for the Floquet equation (|^ . 

We assume a driving that leaves all off-diagonal matrix 
elements of the wire Hamiltonian time-independent while 
the tight-binding levels undergo a position-dependent, 
time-periodic driving fnit) = fn{t + T) with zero time- 
average. Then, the wire Hamiltonian is of the form 

Hwirc(i) = rCo + Y^fnit) \n){n\. (82) 

n 

If HQ, represents the largest energy scale of the prob- 
lem, we can in the Floquet equation (|55|) treat the static 
part of the Hamiltonian as a perturbation. Correspond- 
ingly, the eigenfunctions of the operator fn{t)\n) {n \ — 
ihd/dt determine the zeroth order Floquet states 

e-*^"(*V)- (83) 
We have defined the function 

^n{t) = j:j dt' frXt')^Fn{t + T), (84) 

which is T-periodic due to the zero time-average of fn{t). 
As a consequence of this periodicity, to zeroth order the 
quasienergies are zero (mod ftfi) and the Floquet spec- 
trum is given by multiples of the photon energy, kMl. 
Each fc = 0,±l,±2,... defines a degenerate subspace of 
the extended Hilbert space. If now Ml is larger than all 
other energy scales, the first-order correction to the Flo- 
quet states and the quasienergies can be calculated by 
diagonalizing the perturbation in the subspace defined 
by fc = 0. Thus, we have to solve the time-independent 
eigenvalue equation 

{n.s-ima) = {el-ih-il)\a). (85) 



The time-independent effective Hamiltonian TYcff is de- 
fined by the matrix elements of the original static Hamil- 
tonian TLo with the zeroth order Floquet states (jHSJ, 

(Heff)„„' = P'^"^'\Ho)nn'e~^^-'^'\ (86) 

The i-integration constitutes the inner product in the 
Hilbert space extended by a periodic time coordinatei^ 
To first order in l/fi, the quasienergies — ifvyl^ are 
given by the eigenvalues of the static equation l|S5|l and, 
consequently, the corresponding Floquet states read 

|«„(t)} =^e-'^"(*V>("|a>- (87) 

n 

The fact that all Fn{t) are T-periodic, allows to write in 
(|S7|l the time-dependent phase factor as a Fourier series, 

e-'^"(*) =^a„,fce-^'="*. (88) 
fc 

Thus, {n\ua^k) = o,n.k{n\oi) and the Green function for 
the high-frequency driving reads 

fc' 

where G°^(e) denotes the Green function corresponding 
to the static Hamiltonian Tioff with the self-energy S. 
Finally, substituting e — > e + k'Ml and using the sum 
rule X]fc' an,fc+fc'a* = (5fc,o, we obtain 

I^lj deT,s{e){fRM^)- hM^)]- (90) 

The effective transmission Tcs{t) — V R\G'^^{f)\'^ is 
computed from the effective Hamiltonian H86|l: the elec- 
tron distribution is given by 

/L,eff(e) = ^|ai,fcP/L(e + A:?if^) (91) 
fc 

and /ij,eff follows from the replacement (1, L) — > {N, R). 

In order to derive a high-frequency approximation for 
the zero- frequency noise S, we insert the Green func- 
tion into and neglect products of the type 
C'^ie - kKl) C'^ie - k'Kl) for k ^ k' . Employing the 
above sum rule for the Fourier coefhcients a„^fc, we obtain 
for the noise the static expression Ht)2|) but with the trans- 
mission T(e) and the Fermi functions //{,L(e) replaced by 
the effective transmission Tcff (e) and the effective distri- 
bution function (|91|l . respectively. The fact that /L,cff (e) 
is generally not a mere Heaviside step function has an 
intriguing consequence: In the presence of driving, the 
noise remains finite even if both voltage and tempera- 
ture are zero. 

Two differences between the high-frequency approxi- 
mation and the homogeneous driving, cf. Sec. IIV Dl are 
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worth mentioning: First, the static transmission is now 
replaced by an effective transmission which can be con- 
siderably inffuenced by the driving, see below. Second, in 
general ai^k ^ aN,k such that Jr^cS ^ /l.cS- This means 
that the driving can create an effective bias and thereby 
create a non-adiabatic pump current. By contrast Eq. 
(|80|l reveals that a homogeneous driving cannot create 
such a pump current. Moreover, if all Fn are identical 
as in the case of a homogeneous driving, the effective 
Hamiltonian TLaS equals the original static Hamiltonian. 
Then, also the second line of Eq. (|81|l can be written in 
terms of the static transmission T(e). 

V. CONDUCTOR DRIVEN BY AN 
OSCILLATING DIPOLE FIELD 

In this section, we apply the formalism derived in 
Sees, mil and II VI to study the conduction and noise prop- 
erties of a nanoscale conductor under the influence of 
an electromagnetic field. As an elementary model that 
captures the essential features of a molecular wirCfS. we 
employ a tight-binding model composed of N sites as 
sketched in Fig. ^ Each orbital is coupled to its near- 
est neighbor by a hopping matrix element A, thus, the 
single-particle wire Hamiltonian reads 

JV-l 

n=l (92) 

n 

where En denote the on-site energies of the tight-binding 
levels. Within a dipole approximation, the oscillat- 
ing electromagnetic field causes the time-dependent level 
shifts 

f.n{t) = Acofi{nt)xn (93) 

with Xn = {N + 1 — 2n)/2 the scaled position of site 
\n). Since typical laser frequencies are below the work 
function of a usual metal, we assume that the radiation 
does not penetrate the leads and that, consequently, the 
leads stay in thermal equilibrium. The energy A denotes 
the electrical field amplitude multiplied by the electron 
charge and the distance between two neighboring sites 
and, thus, depends implicitly on the length of the sample. 
This model describes, as well, an array of coherently cou- 
pled quantum dotg^'^''^ under the influence of microwave 
radiation. 

The dipole approximation inherent to the driving H93I) 
neglects the propagation of the electromagnetic field and, 
thus, is valid only for wavelengths that are much larger 
than the size of the samplei^ This condition is indeed ful- 
filled for both applications we have in mind: For molecu- 
lar wires, we consider frequencies up the optical spectral 
range, i.e., wavelengths of the order 1 /im and samples 
that extend over a few nanometers. Coupled quantum 



dots typicallj*ii2i^ have a distance of less than 1 /im while 
the coupling matrix element A is of the order of 30 /leV 
which corresponds to a wavelength of roughly 1 cm. 

We assume that the wire couples equally strong to both 
leads, thus, F^ — Tji = F. An applied transport voltage 
V is mapped to a symmetric shift of the leads' chemical 
potentials, fiu = —^l = e.V/2. Moreover, for the eval- 
uation of the dc current and the zero-frequency noise, 
we restrict ourselves to zero temperature. The zero- 
temperature limit is physically well justified for molec- 
ular wires at room temperature and for quantum dots at 
helium temperature since in both cases thermal electron 
excitations do not play a significant role. 

A. Current and noise suppression 

For a wire described by the Hamiltonian H92() . it has 
been found™ that a dipole force of the form (|93|l sup- 
presses the transport if the ratio A/Ml is close to a zero of 
the Bessel function Jq (i.e., values 2.405.., 5.520.., 8.654.., 
. . . ). Moreover, in the vicinity of such suppressions, the 
shot noise characterized by the Fano factor H13|l assumes 
two characteristic minima. These suppression effects are 
most pronounced in the high-frequency regime, i.e., if the 
energy quanta Ml of the driving exceed the energy scales 
of the wire. Thus, before going into a detailed discussion, 
we start with a qualitative description of the effect based 
on the static approximation for a high-frequency driving 
that has been derived in Sec. IIV El 

Let us consider first the limit of a voltage which is so 
large that in Eq. /i?,cff — fh.cS can be replaced by 

unity. Then, the average current is determined by the 
effective Hamiltonian 

N~l N 
71—1 n—1 

(94) 

which has been derived by inserting the driving (|93|l into 
Eqs. H84|l and (|86|l . Then, obviously HoS is identical 
to the Hamiltonian H92|l in the absence of the driving 
field but with the tunnel matrix element renormalized 
according to 

A ^ Aofi = Jo{A/hn)A. (95) 

Since the Bessel function Jo assumes values between zero 
and one, the amplitude of the driving field allows to 
switch the absolute value of the effective hopping on the 
wire, Acff, between and A. Since the transmission of 
an undriven wire is proportional to |Ap, the effective 
transmission Tcs{e) acquires a factor jQ{A/hn). This 
renormalization of the hopping results finally in a cur- 
rent suppressioni^2i2£ 

For the discussion of the shot noise, we employ the 
Fano factor (|13|l as a measure. In the limit of large 
applied voltages, we have to distinguish two limits: (i) 
weak wire-lead coupling F <C Agff (i.e., weak with re- 
spect to the effective hopping) and (ii) strong wire-lead 
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coupling r ^ Aeff. In the first case, the tunnel con- 
tacts between the lead and the wire act as "bottlenecks" 
for the transport. In that sense they form barriers. Thus 
qualitatively, we face a double barrier situation and, con- 
sequently, expect the shot noise to exhibit a Fano factor 
F « 1/2.'^^ In the second case, the links between the 
wire sites act as — 1 barriers. Correspondingly, the 
Fano factor assumes values F w 1 for N = 2 (single bar- 
rier) and F « 1/2 for TV = 3 (double barrier >S2,. At the 
crossover between the two limits, the conductor is opti- 
mally "barrier free" such that the Fano factor assumes 
its minimum. 

In order to be more quantitative, we evaluate the cur- 
rent and the zero-frequency noise in more detail thereby 
considering a finite voltage. This requires a closer look 
at the effective electron distribution (|91|l : in particular, 
we have to quantify the concept of a "practically infinite" 
voltage. In a static situation, the voltage can be replaced 
by infinity, /^(e) = 1 = 1 — /L(e), if all eigenenergies of 
the wire lie well inside the range [^L,^J'R]■ In contrast 
to the Fermi functions, the effective electron distribution 
(|91|l which is decisive here, decays over a broad range in 
multiple steps of size hfl. Since for our model, Tcs{e) is 
peaked around e = 0, we replace here the effective elec- 
tron distributions by their values for e = 0, 



Jt 



A{N -I] 

2hn 



(96) 



for zero temperature. We have inserted the coefficients 
ai,k = Jk{A{N-l)/2nn) andflAT^fc = J^k{A{N -I) /2nn) 
which have been computed directly from their definition 
(|88|l : Jfc denotes the kih order Bessel functions of the 
first kind. The current, the noise, and the Fano factor 
are given by the static expressions H6U|) and l|62|) with 
the transmission and the electron distribution replaced 
by the corresponding effective quantities, Tcff and /cff.£, 
respectively. Thus, we obtain 



/=A/oo, 



1-A^ 
2A ' 



(97) 
(98) 

(99) 



respectively, where the subscript oo denotes the corre- 
sponding quantities in the infinite voltage limit. 



(100) 



^oo I deT,sie)[l - T,s{e)], (101) 



and Foo — Soo/eloo- The factor 

A = /fl,cff(0) - /L,eff(0) = 

\k\<K{V) 



A{N -I) 

2m 



(102) 



reflects the influence of a finite voltage; KiV) denotes the 
largest integer not exceeding e|y|/2fi,ri. Since Jk{x) « 
for |fc| > X and X^fc'^K^) ~ li find A = 1 if 
K{V) > A{N ~ l)/2hn. This means that for smaU driv- 
ing amplitudes A < eV/{N — 1), we can consider the 
voltage as practically infinite. With an increasing driv- 
ing strength, A decreases and, thus, the current becomes 
smaller by a factor A but still exhibits suppressions. By 
contrast, since Foo < 1 for all situations considered here 
[cf. the remark after Eq. we find from Eq. (O that 

the Fano factor will increase with smaller A. 



B. Numerical results 

The qualitative discussion of the current and noise sup- 
pressions can be corroborated by exact numerical results. 
For this purpose, we have solved numerically the Flo- 
quet equation (|55|) . With the resulting Floquet states 
and quasienergies, we obtained the Green function H58|l . 
In the zero temperature limit considered here, the Fermi 
functions in the expressions for the average current l|37() 
and the zero-frequency noise (|45|) become step functions. 
The remaining energy integrals can be performed analyt- 
ically since the integrands are rational functions. 



1. Intermediate wire-lead coupling 

Figure [21 depicts the average current, the zero- 
frequency noise, and the corresponding Fano factor for 
a wire that consists of iV = 3 sites with on-site energies 
En = as sketched in Fig. ^ The driving frequency 
il — 5A/h lies above all transition energies of the wire 
states and the applied voltage V = 48A/e is relatively 
large. This particular value of the voltage has been se- 
lected to avoid chemical potentials to lie close to multi- 
ples of hfi, i.e., close to the steps of the effective electron 
distribution (|96|l . The wire-lead coupling F — 0.5A is 
sufficiently weak, such that in the absence of the driving, 
the transport is dominated by resonant tunneling. Cor- 
respondingly, the current is essentially determined by the 
hopping rate r/2?i of the electrons from the lead to the 
wire. The noise exhibits a Fano factor F « 1/2 which 
is the characteristic value for the transport across a dou- 
ble barrieri^2iS With an increasing driving amplitude, 
the current becomes smaller until it reaches its minimum 
when the ratio A / assumes a zero of the Bessel func- 
tion Jo- Note that while the analytical treatment within 
a high-frequency approximation predicts a vanishing cur- 
rent, the exact result is still roughly 1% of the value in 
the absence of the driving. Close to the current suppres- 
sion, the effective tunnel matrix element (I95|l is much 
smaller than the wire-lead coupling F and the connec- 
tions to the central site of the wire form a double bar- 
rier. Consequently, we again find a Fano factor F « 1/2. 
At the crossover A^s ~ F, the effective barriers vanish 
and, therefore, the Fano factor assumes its minimum. 
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FIG. 2: Time-averaged current (a), zero-frequency noise (b), 
and Fano factor (c) for a conductor consisting of = 3 sites 
with equal on-site energies, En = 0, as a functions of the 
driving amplitude A. The driving frequency is O = 5A/h, the 
wire-lead coupling is F = 0.5A, and the chemical potentials 
are — —fJ,L = 24A. The exact numerical results (solid 
lines) are compared to the high-frequency approximation for 
finite (dashed) and infinite voltage (dash-dotted). 



These exact numerical results are well reproduced by the 
expressions (|97|l - (|99f) obtained within a high-frequency 
approximations for finite voltage. Figure [3 also reveals 
that for small driving amplitudes, A < eV/{N — 1), the 
assumption of a practically infinite voltage yields the cor- 
rect results. By contrast, for larger driving amplitudes, 
A > eV/{N — 1), the Fano factor can assume values 
F > 1, i.e., the shot noise becomes even larger than in 
the static situation. 

As the coupling strength F is lowered, the distance 
between a pair of minima of the Fano factor becomes 
smaller until the minima finally vanish"^^ (not shown). 
In this limit, the current and the noise are given by the 
weak coupling results (|74|) and H75|l . respectively. The 
corresponding Fano factor = 1/2 [cf. the discussion 
after Eq. I|75|l ] is independent of the driving amplitude. 

Figure |3i depicts the behavior for a driving frequency 
which is of the order of the wire excitations, = A/h. 
Then, the high-frequency approximation is no longer ap- 
plicable. Nevertheless, the average current exhibits clear 
minima with a reduction of the order 50%. Compared 
to the high-frequency case, these minima are shifted to- 
wards smaller driving amplitudes, i.e., they occur for ra- 
tios A/hn slightly below the zeros of the Bessel function 



FIG. 3: (a) Time-averaged current (solid line) and zero- 
frequency noise (dashed line) as a function of the driving 
amplitude for the driving frequency = A/h. (b) Corre- 
sponding Fano factor for the same data (solid line) and for 
the driving frequencies = 1.5A/ft (broken) and = 3A/h 
(dash-dotted). All other parameters are as in Fig.|21 

Jq. At the minima of the current, the Fano factor (solid 
line in Fig. Eh) still assumes a maximum with a value 
close to F w 1/2. Although, the sharp minima close to 
the current suppressions have vanished, in-between the 
maxima the Fano factor assumes remarkably low val- 
ues {F w 0.2). Figure |3Jd also reveals that already for 
51 w 3 A / h, the high-frequency regime is reached. 



2. Strong wire-lead coupling 

For strong wire-lead couphng, it is possible to choose a 
driving frequency that is large with respect to the wire ex- 
citations, but small as compared to the coupling F, thus 
A ^ hft ^ F. Figure0]depicts the current and the Fano 
factor in this limit for wires with a different number of 
sites. The qualitative difference between these cases can 
be explained by the fact that due to the strong coupling, 
the first and the last wire site hybridize with the leads. 
Then the setup behaves similar to a wire with N — 2 
sites and a weak wire-lead coupling cx A^/F. This means 
that for = 2 the wire acts as point contact while for 
A^ = 3, we qualitatively have resonant transport through 
a single level. In both cases remains no tunneling matrix 
element of the wire that could be renormalized and, con- 
sequently, for A^ < 3 the current suppressions vanish in 
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FIG. 4: Time-averaged current (a) and Fano factor (b) as a 
function of the driving amplitude A for a wire with A'' = 2,3,4 
sites and the wire-lead coupling strength V = 10 A. The other 
parameters are En = 0, f2 = bl^/h and = — /ii = 25A. 



FIG. 5: Time-averaged current (a) and Fano factor (b) as a 
function of the driving amplitude A for a wire with = 3 
sites in the presence of an internal bias. The on-site energies 
are Ei — b, E2 = 0, £3 — —b. All other parameters are as in 
Fig.H 



the strong-coupling limit (cf. Fig. This scenario is 
also reflected in the behavior of the Fano factor (Fig. ^p) 
which exhibits the characteristic values F 1 (point con- 
tact) for N = 2 and F « 1/2 (single resonant level) for 
N = 3. Finally, for A^ = 4 we observe the behavior of 
a driven wire with two sites and weak couplingi^ Then, 
a vanishing effective hopping Aoff ~ corresponds to a 
point contact, thus, F Ri 1. Although the behavior of the 
Fano factor can be explained by drawing analogies to a 
weakly coupled wire with N — 2 sites, the global decay 
of the current with the driving amplitude, cf. Fig. is 
not within the scope of this intuitive picture. 



3. Internal bias 

So far, we have assumed that all on-site energies of the 
wire are identical. In an experimental setup, however, the 
applied transport voltage acts also a static dipole force 
which rearranges the charge distribution in the conductor 
and thereby causes an internal potential profilei^2i22i2i 
The self-consistent treatment of such effects is, in par- 
ticular in the time-dependent case, rather ambitious and 
beyond the scope of this work. Thus, here we only de- 
rive the consequences of a static bias without determining 
its shape from microscopic considerations. We assume a 
position-dependent static shift of the on-site energies by 



an energy — &x„, i.e., for a wire with A^ = 3 sites. 

El =6, = 0, £^3 = -b. (103) 

Figure demonstrates that the behavior of the aver- 
age current is fairly stable against the bias. In particu- 
lar, we still find pronounced current suppressions. Note 
that since b ^ ft a high-frequency approximation is still 
possible. As a main effect of the bias, we find reduced 
current maxima while the minima remain. By contrast, 
the minima of the Fano factor (Fig. Et) become washed 
out: Once the bias becomes of the order of the wire- 
lead coupling, 6 « F, the structure in the Fano factor 
vanishes and we find _F w 1/2 for all driving ampli- 
tudes A < eV/ {N — 1) [cf. the discussion discussion after 
Eq. p02|l ] . Interestingly, the value of the Fano factor at 
current suppressions is bias independent. 



VI. CONCLUSIONS 

We have derived with Eqs. and ^ 

expressions for the dc current, the zero-frequency noise, 
and the time-dependent current for the electron trans- 
port through ac-driven nanoscale systems. A corner- 
stone of our approach is the relation of the propaga- 
tor to a non-Hermitian Floquet equation. This yields 
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explicit formulae for the current and the noise. More- 
over, the connection to Floquet theory allows to eluci- 
date various approximation schemes that enable an effi- 
cient computation and, in addition, provide physical in- 
sight. Above all, a high-frequency approximation has 
emerged to be very useful: Within an expansion in 1/0, 
the driven transport problem can be approximated by 
a time-independent transport problem with a renormal- 
ized tunneling and effective distribution functions for the 
lead electrons. The conductance properties of the lat- 
ter can be derived with standard methods. Moreover, 
for the case of a time-dependent gating voltage, we have 
revealed the limitations of the Tien-Gordon approach: 
While such a treatment provides the correct expression 
for the current, it neglects interferences of different Flo- 
quet channels. 

A detailed investigation of the recently found shot 
noise suppression provided a deeper understanding of this 
effect. In particular, the analytical treatment within a 
high-frequency approximation can explain the character- 
istic emergence of the current suppressions which are ac- 
companied by a noise maximum and two remarkably low 
minima. A numerical study fully confirmed the analyt- 
ical results. For lower driving frequencies, i.e., beyond 
the high-frequency limit, the current suppressions be- 
come considerably less pronounced. By contrast, the shot 
noise suppression turned out to be more stable. Thus, 
since the current stays remarkably large while the noise 
is controllable, this regime is particularly promising for 
applications. At first sight, in the limit of strong wire- 
lead coupling these phenomena appear quite different. A 
closer look, however, revealed that the strong coupling 
entails a hybridization of the first and the last site with 
the respective lead. Therefore, the wire behaves quali- 
tatively like a weakly coupled wire with two sites less. 
Moreover, we have found that the noise suppressions are 
quite sensitive to an internal bias. Once the on-site en- 
ergies of neighboring sites have differences of the order 
of the wire-lead coupling energy, the minima of the Fano 
factor vanish. A most interesting application of these 
results is the development of current sources with a con- 
trollable noise level. 
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APPENDIX A: AC TRANSPORT VOLTAGE 

Within this work, we focus on models where the driv- 
ing enters solely by means of time-dependent matrix 
elements of the wire Hamiltonian while the leads and 
the wire- lead couplings remain time- independent. An a 
priori different type of driving is the application of a 
time-dependent transport voltage. In this appendix, we 
demonstrate that a setup with an oscillating transport 
voltage can be mapped by a gauge transformation to a 
Hamiltonian of the form . Consequently, it is possible 
to apply the formalism derived derived in Sec. IIIII also to 
situations with an oscillating transport voltage. 

We restrict the discussion to the case where the elec- 
tron energies of only the left lead are modified by an ex- 
ternal T-periodic voltage Vac(i) with zero time-average, 
thus in the left lead 

eq^eg-eV^t). (Al) 

The generalization to a situation where also the lev- 
els in the right lead are T-periodically time-dependent, 
is straightforward. Since an externally applied voltage 
causes a potential drop along the wire)S2i22i2i we have to 
assume for consistency that for an ac voltage, the wire 
Hamiltonian also obeys a time-dependence. Ignoring 
such a time-dependent potential profile enables a treat- 
ment of the transport problem within the approach of 
Refs. 63 and 64. In the general case, however, we have 
to resort to the approach put forward with this work. 

We start out by a gauge transformation of the Hamil- 
tonian with the unitary operator 

[/ac(t) = exp { - ic^it) (clci + J2 Aq^L^ } (A2) 

9 

where 




describes the phase accumulated from the oscillating 
voltage. The transformation ljA2(l has been constructed 
such that the new Hamiltonian H[t) = Ul.H{t)U^c - 
iKUl^Jj^^ possesses a time-independent tunnel coupling. 
Since, the operator ci transforms as ci ci exp(— 2(/)(t)), 
the matrix elements -ff„„' {t) of the wire Hamiltonian ac- 
quire an additional time-dependence. 

The second term in the Hamiltonian ljA4|l stems from 
—ihUl^U^^. Owing to the zero time-average of the volt- 
age Vac(i), the phase (f){t) is T-periodic. Therefore, 
the transformed wire Hamiltonian is also T-periodic 
while the contact and the lead contributions are time- 
independent, thus, H{t) is of form 
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APPENDIX B: ALTERNATIVE DERIVATION 

In Ref-ls^, the expressions and for the current 
and the noise in the wide-band hmit have been derived 
by eHminating the leads in favor of a stochastic opera- 
tor. In this appendix, we detail this approach. Like in 
Section irm we start here also from the Heisenberg equa- 
tions H14(l - H16(l for the annihilation operators. The ones 
for the lead operators, Eq. 1)14(1 . are easily integrated to 
read 

t-to 


(Bl) 

and C}iq{t) accordingly. Inserting (|Bip into the Heisen- 
berg equations 1(15(1 for the wire operators yields 



Cl/N 



T-y^^l/Af.ri'ft)' 



L/R 

~2h' 



Cl/N + iL/R{t), 



c„ = --^^i/„„'(i)c„', n = 2,...,7V-l. (B2) 

n' 

Owing to the wide-band limit, the dissipative terms 
are memory free. Within the chosen grand canonical 
ensembles the operator- valued Gaussian noise S,e,{t) — 
-(*/^) E„ "t^;„e— «(*-*°)/''Q,(to) obeys 



m)) = 0, 



(B3) 



(C|(t)6'(t')) = /dee"(*-*')/''/,(e) . (B4) 



The current operator then assumes the form 
Idt) = |ricl(i)ci(i)-e{cl(t)eLW+dWciW}- (B5) 

The homogeneous set of equations that corresponds 
to ((B2ll coincides with the equations of motion ((46() and 
((50(1 which are solved by the Floquet states \ua{t)). Thus, 
the Floquet states \ua{t)) together with the adjoint states 



\u^{t)), allow to write the solution of l(B2(l in closed form. 
In the asymptotic limit Iq — > — oo, it reads 

Cn{t)= dT{n\U{t,t-T) 

Jo (B6) 

X {rndt ~ r) + \N)Ut - r)} . 

where U{t,t — r) is the propagator l(56() for the wire elec- 
trons. 

To obtain the current {Ii^{t)), we insert the operator 
((B6(l into the expression ((40(1 and use the expectation 
values ((B4ll . With the Green function l(S^ . we find the 
still unsymmetric expression 

lde{rL\Gn{t,e)\'fL{e) 
+TR\GiN{t,e)\^fR{e) 
+i[Gl,{t,e)-Gn{t,e)]fL{e)}. 

For a symmetrization, we eliminate the backscattering 
terms, i.e., terms containing Gn, by use of the relation^ 



G\t, e) - G{t, e) = ^h-G\t, e)G{t, e) ^^^^ 
+ 2iG\t,e)j:G{t,e) 

which follows readily from the Floquet representation 
((56(1 of the propagator and the Floquet eigenvalue equa- 
tion l(55() together with its adjoint. A subsequent Fourier 
transformation with respect to t — t~t' yields Eq. ((B8(l . 
By inserting the matrix element (1 1 . . . 1 1) , we obtain from 
((B7() for the time-dependent current the symmetric ex- 
pression l(4U() . 

To derive an expression for the zero frequency noise, 
we insert the operator l(B6() into the definition (|5Jl of the 
current-current correlation function and integrate over 
the times r and t. Again, we employ the relation l(B8() to 
bring S into the symmetric form ((45(1 . 
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